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Abstract 

We show that in the PQCD approach the B pOK decays display not only the dynam¬ 
ical penguin enhancement, the mechanism responsible for the large branching ratios (BRs) 
of -B —> cj)K decays, but also the importance of annihilation contributions. We hnd that the 
CP asymmetries (CPAs) of 5 —> p'^K^, —>• p^{uj)K^ (class I) are all over 50%. 


PACS number(s): 12.38.BX, 12.38.QK, 13.25.HW 


The study of charmless B meson decays has an enormous progress, since many modes, 
such as B ^ 7171 and B tiK, were measured by CLEO |I[, as well as by BABAR [Q and 
BELLE 0. Via the search of B decays, we can not only test the origin of CP violation in 
standard model (SM), which is the consequence of the Cabibbo-Kobayashi-Maskawa (CKM 
) quark-mixing matrix Q, but also verify various QCD approaches, proposed to deal with 
hadronic effects in exclusive decays. Recently, the BRs of the modes involving a vector (V) 
meson, such as B ^ measured by CLEO 0 and B^ —> cjyK^ by the B factories, 

B^ (t)K^ = ( 5 . 5 ;?;^ ± 0.6) X 10“® (CLEO §), 

= (7.7;};®±0.8) X 10“® (BABAR 0), 

= (11.2l|2 ± 1.4) X 10"® (BELLE g) 


have been also observed. Although the present results of B^ —»• (j)K^ among the B factories 
do not match well each other, more accurate measurements will be obtained in the near 
future. Nevertheless, the theoretical predictions are not consistent either; for instance, the 
modihed perturbative QCD (MPQCD) predicts (10.2^2;?) x 10“® g, |^, but QCD factor¬ 


ization (QCDF) [0 gives x 10’ -» H, Therefore, an interesting question is raised: 

besides B —>■ (fiK decays, can we hnd other processes, that distinguish which QCD approach 
is proper and correct? or what processes can support either of them? Following the spec¬ 
ulation oi B ^ VP decays, P being a pseudoscalar meson, we hnd that it is important to 
investigate B V'K decays {V' = p and a;), because many properties, such as CPAs of 
over 50% and essential annihilation contributions, are remarkable. Although the V meson 
carries spin degrees of freedom, only the longitudinal component contributes. Thus, in order 
to understand the properties of R —> V'K, it is useful to combine the analysis oi B ^ 7r~K^ 
and B^ —> decays, all of which are described by the b —>■ s transition. Note that the 

pion can be regarded as the longitudinal component of a p meson. 

The effective Hamiltonian for decays with the 6 —>■ s transition is given by 


ft. = % E n 

V ^ q'=u,c 


10 

+ c'2(/i)(P?'^ + J2c^ip)o, 


1=3 


( 1 ) 


where Vq' = V*,gVqib are the products of the CKM matrix elements, Cj(/i) are the Wilson 
coefficients (WCs) and Oi correspond to the four-quark operators. The explicit expressions 
of Ci{p) and Oi can be found in Ref.[^. We can dehne the evolution variables as 
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where the superscripts q represent the light u, d and s quarks. Note that 02 is larger than 
oi and nonfactorizable effects are only associated with the color suppressed variables 

For a simple analysis, we hrst concentrate on R — p~ decay, because it contains the 
common properties of R —> V'K decays. Hence, based on Eq. (|^), the decay amplitudes for 
R —> 7i~K^, R+ —»• (pK^ and R —>■ p~K^ are written as 

v; (al“^ + 2r^4“)) - K*«2 
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with tk = mPj^/MB, rnPx being the chiral symmetry breaking (CSB) parameter associated 
with the kaon. We present only the dominant and subdominant effects with nonfactorizable 
effects neglected. However, the complete effects will be included in our numerical calcula¬ 
tions. fx, f(i> and fs are the decay constants of K, 0 and B mesons, respectively. 
denote the B —>■ P{V) form factors, and are the time-like form factors induced by 

the annihilation contributions from the (S' — P)(S' -I- P) four-quark operators. As for 
generated by the {V — A){y — A) operators, due to helicity suppression, we do not discuss 
them. We emphasis that the expressions of in Eqs. (^-(^ are given only for 

the short-hand convenience. In the MPQCD approach, they are convolutions of and 
^e[a] ^ actually. If we apply the idea of QCDF to MPQCD, in the viewpoint of numerical 
estimation, the /i scale of could be hxed at around the scale of 1.7 GeV. The scale is 


around 2.7 GeV in QGDF 

Although C 2 is one order of magnitude larger than G 4 and Gg, which are the hrst two 
largest WGs of QGD penguin, the penguin emission topologies in the tt~K^ and modes 
are the dominant effects due to the suppression of GKM matrix element W ~ 1.6AA®e“*'^. 
The tree and penguin annihilation are subdominant. If taking 7 ~ 90°, the numerical value 
of Eq. d^) is larger than that of Eq. (|) by an approximate factor of a/L8 and a/^ with /i 
being fixed at the different scales in MPQGD and QGDF, respectively. As known, the world 
average BR of P —> is (17.2 ± 1.5) x 10“° JI], ||, |^. Following the above analysis, we 

expect the BR of (pK^ mode is around 9.5 x 10“° (/i = 1.7 GeV) and 5.5 x 10“° (/i = 2.7 
GeV). Glearly, with the information of P —in the present experimental status, we 
cannot determine the typical QGD scale for heavy-to-light decays, unless more precision 
measurements on P —>■ (pK are given. 

Different from P —> 7 r“A'’'' and B^ —^ (pK~^ decays, in terms of Eq. (|^), we see that due to 
the cancellation between the two terms in — 2rKa'i\ the penguin emission contributions 
of P —>■ p~K^ decay are no longer leading. Instead, tree diagram may be the main effects. 
It is known that a 2 {p) is a flat function in p. However, af^{p) and a^\p) have an enormous 
enhancement aX p < Mb/2 [M. If F^^ is only few times smaller than Pf^, with the 
penguin enhancement from (p), the annihilation contributions could possibly become 
the essential parts in P ^ p~K^ decay. Recently, the estimation of annihilation topologies 
has been proposed by MPQGD and QGDF |^. Although the basic idea is originated 


from the Lepage and Brodsky (LB) formalism |Q, in which a transition amplitude can be 


factorized into the convolution of nonperturbative parts, described by hadron wave functions, 
and a hard amplitude of valence quarks, dictated by perturbative hard gluon exchanges, both 
approaches have different description on factorization scale. 

It has been known that the original LB formalism suffers singularities from the end¬ 
point region with a momentum fraction a; —> 0. In the QGDF approach, it is claimed that 
heavy P meson decays at the fast recoil region are still dominated by soft gluons, which 
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are uncalculable perturbatively. As a result, the decay amplitudes for B meson decays 
are written as A ~ with to ~ Mb ~ Mb/2, where C{to), /mi and 

are the relevant wilson coefficient, decay constant of Mi meson and B M 2 form factor 
respectively. However, following the concept of PQCD, if the spectator quark of inside B 
meson, carrying the momentum of order of A with A = Mb — M^, and M^ being b quark 
mass, wants to catch up with the outgoing quark, which is the daughter of b quark decay, 
to form a hadron, it should need to obtain a large energy from b quark or the daughter 
of b quark. That is, in contrast to the conclusion of QCF approach, hard gluons actually 
play an essential role in B meson decays. Therefore, the relevant decay amplitude should 
be calculable perturbatively and the order of magnitude of typical scale is similar to the 


momentum of hard gluon approximately, denoted by y AMb 

Now, how to deal with the problem of singularities is the main part of PQCD. In order to 
handle these singularities, the strategy of including kT, the transverse momentum of the va¬ 
lence quark, and threshold resummation has been proposed . It has been shown that the 
singularities do not exist in a self-consistent MPQCD analysis [ 0 - Additionally, MPQCD 
gives many interesting results: for example, the B P{V) form factors are dominated by 
perturbative dynamics with Os/tt < 0 .2, the formalism involves less theoretical uncertain 
parameters that arise from the shape parameter ujb in the B meson wave function, CSB 
parameter mPj^ and the parametrization of Sudakov factor from threshold resummation, the 
penguin contribution is enhanced at a lower typical scale and annihilation topologies con¬ 
tribute large absorptive parts. Especially, according to our power counting rules, the ratios 
of the transition form factor to the annihilation contributions and to the nonfactoriz- 

able contributions are found to be = 1 : A/Mb m 


For Mb ~ 5.0 GeV, the magnitude 
of the annihilation amplitude is only less than that of F^^ by a factor of 3. In the literature, 
the applications of the MPQCD to the processes of B 

K'qb) 


g, H ^ KK 


B 


and Ho —KK 


PP, such as B ^ Kn |]^, B 
, that oi B VP such as B 


TTTT 

(/n 


^ and B |Q, as well as that of semileptonic decays such as H — K^*H~^i~, have 

been studied and found that they are consistent with the experimental data or limits. 
Hence, the H —>■ p and time-like form factors, dehned as 


< P {PpPp)\b AlbulBipB) > = 2mpe*-qFfA 

{p~{Pp)K^{Pk)\uj5s\ 0) = 2mpe*p-qFf 

with q = Pb — Pp, can be explicitly written as 

p 1 poo 
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( 0 p {X 2 ) + (/>p (xs))] He K {xi ,X2,bi,b2) 
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2rp4>^p{x2) He ( 4 ^^) K{x2,Xi,b2,bi)"^ 


Xa6 


= —S'kCfM% / dx2dx3 / 62 ^ 6253^63 

Jo Jo 

{[xsTkc/p {X2) {(/k {xs) - c/k (xo)) 

+ 2 rp(/^p (X2) c/k (a^s) H„ ha {x2, X3, 62, 63) 

+ [-rpX2(/K (xs) (0* (xa) - c/p (xa)) 

+ 2 rK(j)p (X2) c/k (a^s)] H^ ha {x^, X2, 63 , 62 )} 


( 6 ) 


( 7 ) 
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where Vp = mp/Ms, rap is the p meson mass, Cp = 4/3 is the color factor, Xi {i = 1, 2, 3) are 
the momentum fractions carried by the spectator quarks inside the B, p and K mesons, and 
the variables bi are conjugate to parton transverse momenta kip- (f)p(v){x) denote the twist- 
2 P(V) meson wave functions, and 0 P(v)(a;) and ‘i^p(v)(^) correspond to the pseudoscalar 
(tensor) and pseudotensor (scalar) twist-3 wave functions, respectively |^, The wave 
functions are pure nonperturbative effects and universal. In our numerical calculations, we 
will adopt the results of Refs. |^, ^ derived from QCD sum rules. The explicit expressions 
of hard functions he(a) are shown as 


he{Xi,X 2 ,bi,b 2 ) 


= Ko {^yXlX 2 MBbl) St{x 2 ) 

X [e{bi - b2)Ko {^/^MBbi) lo (^^^5^2) 
+0{b2-bl)Ko{^/^MBb2)Ioi^/^MBbl)] , 


( 8 ) 


ha{x 2 ,X 3 ,b 2 ,b 3 ) = 


X 


{sJX2X3MBb2) Stixs) 

0(62 — &3)P^0 ^ (a/ 5 ^^S& 2 ) Jq (v^^S^s) 


+ 6 {b 3 — & 2 )-f^o ^ Jo (\/^^B& 2 ) 


(9) 


where St{x) is the evolution function from threshold resummation parametrized by St 
Nt [x(l - x)Y 


given by 


T5|| , and Kq, Jq, and Jq are the Bessel functions. The evolution factors are 
Ee(a){t) = ttsit) exp[-SB{K)it, Xp3)) - Sp{t,X2)] , (10) 


we obtain Ffp = 0.371^1 

ipB-K _ n q+0.03 
Pg — U-'J- 0.02 


which contain the Sudakov factor from kp resummation ® 0 - 

Adopting fB = 0.19, fp = 0.20, fp = 0.16 GeV and with the allowed values of parameters. 

Via the same procedure, we get = 0 . 35 lo;o 2 and 

All of these form factors are consistent with those from light-cone 
QCD sum rules (LCSR) and quark model (QM) [^. Similarly, the time-like form factor 
is given by F^^ = (—0.46 -|- 7.4i) x 10“^, in which the large absorptive part arises from the 
on-shell condition of internal light quarks. As expected, the magnitude of F^p is only larger 
than that of F^^ by a factor of 5 in the MPQCD approach. In terms of the above results, 
from Eq. (^, if choosing 7 ~ 90°, the contributions of tree and penguin annihilation are 
comparable but are opposite in sign, such that the BR of R ^ p ~is small. However, via 
Vu instead of V*, we hnd immediately that the CP conjugate mode B p^K~ is enhanced. 
It is interesting to question how large the BR can be pushed up by such enhancement. For 
estimation, we assume F^^ ^ F^^ and use the above obtained values, from Eqs. (^) and 
(I), the relation |A(R+ —>• 0iP^)| ~ 1.17 \A{B p'^K~)\ for 7 ~ 90° and for p ^ 1.7 
GeV is obtained. With BR{B^ —»• Ri 10 x 10“°, we get BR{B —> p'^K~) ~ 7 x 10“°. 

Furthermore, due to the destruction between the tree and penguin amplitudes in R —> p~K^ 
decay, one also expects the CPA in R —> p'^K^ is large. 

We have studied the properties of R —>• p^K^ and found that by the constructive inter¬ 
ference between the tree and penguin annihilation topologies, the BR and CPA estimated 
in the MPQCD are not as small as those expected in Ref. |^. As mentioned before, 
basically all the modes in R —> p{uj)K have the similar properties. In addition, for those 
decays involving p° and oj mesons that are composed by {uu — dd) /and {uu F dd) /a/2 , 


respectively, the new terms 


,0 

^35 


■gy.®* and 


,0 

^35 


+ with = 


.( 1 ) 


+ 


will be induced. Because the values of 03 ^^ and are much smaller than those of the 
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properties discussed before will not change. Furthermore, we can separate these decays into 
two classes. Class I consists oi B ^ ^ and —>• ooK^, while B 

B —>■ ujK^ and B^ —>■ are in class II. We hnd that the tree contributions in class I are 

related to a 2 F^P, but in class II they are associated with aiF^^ or a 2 F^^. As known, oi 
is one order of magnitude less than 02 and F^^ is much smaller than F^^q , so that the tree 
effects associated with V*aiF^^ and V*a 2 Fl^^ are quite small. Thus, we expect that the 
large CPAs and BRs exist only in the modes of class I. We emphasize that the dynamical 
variables in the MPQCD approach have to be as a evolution factor convoluted with 

Fe\a^'^^■ formulas can be easily obtained by inserting the relevant a!f\ijL) to Eq. 0 ) 

and setting ^ = t. In our following numerical calculations, we will adopt such formulas. 

With the data of I4s — Vts ^ —AA^, Vub ~ A ^ 0.81, A ~ 0.22, Rb ^ 0.36, 

Nt ~ 1.77, c ~ 0.3, = 1.7 and ujb = 0.4 GeV as the input values and considering all 

factorizable and nonfactorizable effects, the results of BRs with CP average and CPAs are 
shown in Table |^. Note that the CPAs in class I are large. That means the BR of one of 
decay modes much larger than that of its CP conjugate mode. For illustration, they are 
given as {B —> p'^K~, B~ —>• p^K~, B~ —>• ujK~) = {8.72, 3.81, 5.10) xl0“® for 7 ~ 72°. 
We also plot the BRs as a function of 7 in Figure |l|. It is worth of pointing out that 
except tree effects, the B^ and B p~K^ decays have the similar contributions 

from the penguin topologies. According to our estimation, the effect of penguin emission 
associated with — 2rKa^^'^ F^p on the BR is only around 19%. Therefore, the value 
of BR{B^ —> p^K^) in Table |I] almost indicates the effects of annihilation contributions. 
In MPQCD, the uncertainties of hadronic effects can be from the power factor c for the 
parametrization of threshold resummation, ujb and With the allowed regions []T|, we 
hnd that the uncertainties of power factor c and ujb are below 10 % and the error of m%- 
is 20%. The most challenge contributions are from higher order corrections. In order to 
estimate their size, we hx all free parameters to specihc values, shown in the beginning of 
this paragraph, and then change the chosen condition of t scale by a deviation of 30%. We 
hnd that the inhuence on BRs is also around 30%. 

We now give a brief discussion on the nonfactorizable ehects. It is known that besides 
contributions from 04 to the nonfactorizable diagrams of the hard gluon exchanges between 
the valence quarks in p and B mesons, those from 035 ^ — 035 ^ and 03 ^^ + 035 ^ for p° and oj 
modes will also be introduced. According to Eq. (^), the former is only associated with the 
small WCs of Ct^io{p) but the latter is related to C' 4 , 6 . Thus, we expect that the inhuence of 
nonfactorizable ehects on ujK modes is extraordinary. From our estimations, by neglecting 
all nonfactorizable contributions, we hnd that the Brs of R —>■ a;77^(77°) are reduced about 
50%, while the inhuence on Br^pK"^) and Br{pK^) are around 35% and 18%, respectively. 

Finally, we give a remark on the QCDF approach. According to the analysis of Ref. [|1^, 
the BR of B^ (j)K^ can reach 7.3 x 10“°, if the included annihilation ehects are almost 
real. If taking this value as the common result of MPQCD and QCDF, with the similar 
procedure, the values of the CP-averaged BR for B p{yj)K in QCDF can be estimated 
as (R ^ R^ ^ p°77=F, R^ ^ ujK^)^ (2.51, 0.62, 1.28) xlO”® and (R^ ^ p^K^, 

B —> p^K^, B —>■ a;77°)~ (1.63, 1.99, 0.64) xl0“° for 7 ~ 72°. The results in MPQCD are 
35% smaller than those in Table |^. Though some results are similar in both approaches, the 
CPAs from QCDF for class I are found to be only at the few percent level. If the considered 
annihilation ehects are almost imaginary, the CPAs from QCDF for class I can be as large as 
those from MPQCD. However, the BRs from QCDF become much smaller than those from 
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MPQCD. Altogether, even if MPQCD and QCDF have the common prediction for B —> (j)K, 
it is clear that B p{oj)K can distinguish the different QCD approaches. 

In summary, we have performed the analysis of the B —> p{uj)K decays. Although the 
unique dynamical penguin enhancement in MPQCD can be verified by the measurements of 
B (j)K, which distinguish MPQCD from QCDF, B p{oj)K can display not only such an 
enhancement but also the importance of annihilation topologies, especially their absorptive 
parts. We also show that nonfactorizable effects play an important role in i? —> ujK decays. 
Moreover, due to the cancellation in 04 '^^ — 2rK(ie for penguin emission topologies, any siz¬ 
able new physics, with or without new weak CP violating phases contributing to the penguin 
operators O 3 ~ (Pio, will have a remarkable influence on the BRs and CPAs of the modes in 
class I and II. B —>• p{uj)K are also sensitive to physics beyond standard model. Recently, 
CLEO has reported results of R — p(uj)K, such as BR{B —>• p'^K^) = {16.01lIa ± 2.8)xl0“® 
and BR{B^ —> ujK^) = {3.2t.i\g ± 0.8)xl0“® [Q. It is obvious that the central value of 
BR{B^ —>■ ujK^) is the same as our prediction. The further confirmation will be made, 
when the data from the B factories are announced in the near future. 
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Figure 1: BRs of the B p{uj)K decays as a function of the angle 7 . The solid, dashed 
and dash-dotted lines are for (a) Bd B^ P^K^-, and and for (b) 

Bd p^K^, B^ — p^K^, and Bd uK^, respectively. 


Table 1: BRs (in units of 10 ®) with CP average for class I and II respectively and CP 
asymmetries (%) in the B p{oj)K decays for 7 ~ 72°. Nonfactorizable effects are included. 


Class I 

Mode 

B - 


B^ p^K^ 

B^ 

ojK^ 

BR 


5.42 

2.18 


3.22 

CPA 

- 

-60.76 

-74.95 

- 

-58.21 

Class II 

Mode 



B p^K^ 

B - 


BR 


2.96 

2.49 


2.07 

CPA 


-3.74 

-6.18 


6.39 
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